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I.  INTRODUCTION 


Under  our  previous  contract  (N00019-77-C-0299)  an  aircraft  simu¬ 
lation  model  was  developed  which  treated  the  fuselage  as  a  composite 
elliptic  cylinder  and  the  wings  and  stablizers  as  finite  flat  plates. 
That  numerical  solution  was  based  on  a  near  field  formulation  as  pre¬ 
sented  in  Reference  1  such  that  one  could  compute  the  pattern  at  a 
finite  range.  As  the  previous  solution  was  used  to  analyze  various 
private,  commercial,  and  military  aircraft,  it  became  apparent  that 
it  had  several  shortcomings  which  had  to  be  corrected  before  it  would 
be  really  useful  for  a  general  aircraft  problem.  The  main  objection 
with  the  previous  solution  was  in  its  inability  to  simulate  arbitrary 
structures  which  might  be  attached  to  the  basic  aircraft.  For  private 
aircraft,  there  was  no  way  with  the  previous  model  to  simulate  wing 
tanks,  engine  housings,  propellors,  etc.  For  commercial  aircraft,  one 
could  not  treat  "T"-tail  aircraft,  jet  engine  housings,  etc.  For  mili¬ 
tary  aircraft,  one  can  easily  visualize  numerous  structures  being  at¬ 
tached  to  the  basic  fuselage  and  wing.  For  these  reasons,  the  previous 
numerical  solution  of  Reference  1  is  to  be  extended  under  the  present 
contract  to  treat  multiple  plates  which  can  be  attached  to  the  fuselage 
or  not,  as  well  as  being  attached  or  not  to  other  plates.  Using  this 
multiple  plate  approach,  one  can  treat  wing  tanks,  engine  housing,  T- 
tails,  etc.  simply  by  boxing-out  the  structure  using  finite  flat  plates. 
The  finite  plate  models  were  chosen  to  simulate  these  structures  in 
that  they  are  easy  to  input  in  the  code  and  are  much  more  efficient 
to  analyze.  The  validity  of  the  flat  plate  simulations  will  be  studied 
under  the  present  contract  as  well  as  others  which  might  wish  to  use  the 
code  to  treat  a  special  application.  Note  that  under  this  contract, 
the  general  near  field  aircraft  code  is  being  developed  and  tested  for 
simple  models.  Specific  applications  of  the  code  will  be  done  as  the 
need  arises  on  various  other  contracts.  For  example,  it  has  been  used 
to  analyse  private  aircraft  configurations  for  NASA  as  reported  in  Ref¬ 
erence  2. 
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II.  VARIOUS  GTD  TERMS 


In  order  to  illustrate  how  the  present  solution  is  developed  based 
on  our  analysis,  a  set  of  problems  are  presented  which  illustrate  the 
different  GTD  mechanisms  included  in  the  solution.  The  first  example 
illustrates  the  GTD  mechanisms  used  in  our  previous  analysis.  For  this 
purpose  a  single  flat  plate  is  attached  to  a  circular  cylinder  as  shown 
in  Figure  1.  The  various  terms  treated  using  the  previous  solution 


Figure  1.  A  flat  plate  attached  to  a  circular  cylinder. 

are  shown  in  Figure  2.  Note  that  each  pattern  is  normalized  to  the 
same  level  such  that  one  can  see  the  relative  significance  of  each  term. 
In  order  to  better  observe  how  the  total  Geometrical  Theory  of  Diffrac¬ 
tion  (GTD)  solution  evolves,  various  terms  are  superimposed  in  Figure 
3.  Note  that  the  double  diffraction  term  is  being  introduced  to  the 
solution  under  the  present  effort.  An  interesting  result  is  shown  in 
Figure  3a  where  the  source  and  reflected  fields  are  superimposed.  These 
two  terms  form  the  classical  "Geometrical  Optics"  (GO)  solution.  How¬ 
ever,  one  should  note  that  the  GO  solution  is  far  from  complete  as  can 
be  observed  from  the  discontinuities  in  the  pattern.  Our  final  GTD 
solution  shown  in  Figure  3d  is  seen  to  be  continuous  and  should  agree 
with  experimental  as  verified  by  several  examples  given  in  Reference 
1. 
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(c)  diffracted  field  (0) 


(d)  S  +  R  +  D 


Figure  2.  Radiation  patterns  due  to  various  ray  mechanisms. 


(a)  S  +  R 


(b)  double  diffracted 
field  (D/D) 


(c)  S  +  R  +  D 


(d)  S  +  R  +  D  +  D/D 


Figure  3.  Radiation  patterns  due  to  various  combinations 
of  the  GTD  terms. 


The  next  example  is  used  to  illustrate  the  additional  GTD  mech¬ 
anisms  added  under  the  present  contract.  This  geometry  is  shown  in 
Figure  4  in  which  a  second  plate  is  added  to  the  model  discussed  in 


Figure  4.  A  bent  plate  attached  to  a  circular  cylinder. 

the  previous  paragraph.  The  source,  reflected,  and  diffracted  terms 
used  in  our  previous  solution  are  shown  in  Figure  5.  The  superimposed 
pattern  using  these  three  GTD  terms  is  illustrated  in  Figure  5d.  One 
should  note  the  discontinuities  in  this  pattern  at  £  14°  and  45°.  These 
discontinuities  are  compensated  for  by  the  higher-order  GTD  interaction 
terms  illustrated  in  Figure  6.  As  before,  these  terms  are  plotted  rela¬ 
tive  to  the  same  signal  level  such  that  one  can  observe  their  relative 
significance.  It  is  clear  that  these  higher-order  terms  can  be  signifi¬ 
cant  in  certain  sectors  of  the  pattern.  In  order  to  demonstrate  how 
the  complete  GTD  solution  creates  the  desired  total  pattern,  various 
combinations  of  terms  are  illustrated  in  Figure  7.  The  GO  solution  is 
shown  in  Figure  7b  and,  as  before,  several  discontinuities  exist  in 
this  pattern.  The  total  superposition  of  the  recently  introduced  higher- 
order  terms  is  shown  in  Figure  7c.  It  is  very  interesting  that  the 
patterns  illustrated  in  Figure  7b  and  c  superimpose  to  give  the  total 
pattern  shown  in  Figure  7d.  As  with  any  GTD  solution  for  a  complex 
structure,  one  can  compute  higher-  and  higher-order  terms.  However, 


(c)  diffracted  field  (D) 


(d)  S  +  R  +  D 


Figure  5.  Radiation  patterns  due  to  various  GTD  terms. 


(c)  diffracted/ref lected 
field  (D/R) 


(d)  diffracted/diffracted 
field  (D/D) 


Figure  6.  Radiation  patterns  due  to  the  second 
order  interaction  GTD  terms. 


(c)  second  order  interaction 
GTD  terms. 

(R/R  +  R/D  +  D/R  +  D/D) 

Figure  7.  Radiation  patterns  due  to 
of  the  GTD  terms. 


there  are  two  major  problems  with  adding  such  terms  beyond  those  included 
here:  1)  The  accuracy  of  the  GTD  diffraction  solutions  become  ques¬ 
tionable,  2)  The  numerical  solution  becomes  very  inefficient.  Furthermore 
based  on  cases  examined  to  date,  it  appears  that  the  major  structures 
found  on  aircraft  can  be  adequately  solved  using  the  interaction  terms 
considered  here.  In  order  to  illustrate  the  various  GTD  mechanisms 
they  are  all  shown  in  Figure  8  in  terms  of  ray  paths. 

The  previous  example  was  used  to  illustrate  the  various  higher-order 
terms  used  in  the  present  analysis.  This  next  problem  is  used  to  show 
the  various  GTD  terms  for  a  more  realistic  geometry  as  shown  in  Figure 
9.  In  this  case,  additional  plates  are  added  which  might  be  used  to 


Figure  9.  Flat  plate  model  for  the  wing  and  the  engine. 

simulate  the  engine  housing  for  a  private  aircraft.  The  first-order 
terms  (source,  reflected,  and  diffracted)  are  shown  in  Figure  10.  The 
second-order  interaction  terms  are  shown  in  Figure  11,  and  they  repre¬ 
sent  the  modification  to  the  analysis  made  under  the  present  contract. 
Finally,  various  combinations  of  GTD  terms  are  illustrated  in  Figure 
12.  Again,  it  is  very  interesting  how  all  the  different  discontinu¬ 
ous  GTD  terms  combine  to  give  the  nice  smooth  total  pattern  shown  in 
Figure  1 2d . 
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(a)  source  field  (S) 


(b)  reflected  f ield  (R) 


(c)  diffracted  field  (D) 


(d)  S  +  R  +  D 


Figure  10.  Radiation  patterns  due  to  various  GTD  terms. 


(c)  diffracted/reflected  (d)  diffracted/diffracted 

field  (D/R)  field  (D/D) 

Figure  11.  Radiation  patterns  due  to  the  second  order 
iteration  GTD  terms. 


(a)  S  +  R 


(b)  geometric  optics 
solution 
(S  +  R  +  R/R) 


(c)  second  order  interaction  (d)  total  solution 

uiU  terms 

(R/R  +  R/D  +  D/R  +  D/D) 

Figure  12.  Radiation  patterns  due  to  various 
combinations  of  the  GTD  terms. 


III.  EXPERIMENTAL  VERIFICATION 


Since  new  structures  are  being  treated  under  the  present  contract, 
some  experimental  verification  is  necessary.  As  stated  previously, 
this  verification  is  only  done  using  the  simpler  multiple  plate  structures. 
No  attempt  will  be  made  here  by  us  to  measure  scale  models  of  actual 
aircraft  configurations.  This  is  merely  an  attempt  to  verify  that  our 
numerical  solution  is  accurate  for  the  simpler  multiple  plate  structures. 
Comparison  with  scale  model  measurements  will  be  done  on  other  contracts 
such  as  the  private  aircraft  simulations  done  in  Reference  2. 

The  first  example  is  used  to  verify  our  numerical  solution  for 
the  last  configuration  of  the  previous  section.  The  calculated  and 
measured  roll  plane  patterns  are  shown  in  Figure  13.  In  the  next  case, 
the  engine  plate  model  is  slid  out  the  wing.  The  measured  and  calculated 
roll  plane  patterns  for  this  geometry  are  shown  in  Figure  14.  In  each 
case,  excellent  agreement  between  measured  and  calculated  results  is 
obtained. 

Th°  final  example  considered  here  is  illustrated  in  Figure  15. 

This  geometry  is  included  in  that  it  tends  to  resemble  a  more  realistic 
aircraft  configuration.  The  roll  plane  measured  and  calculated  pat¬ 
terns  are  shown  in  Figure  16.  Two  sets  of  measured  and  calculated  pat¬ 
terns  are  illustrated  in  Figures  17  and  18  for  68°  and  112°  conical 

pattern  cuts.  In  each  case  both  the  E  and  E  components  of  the  near 

0  (f) 

zone  field  are  plotted  relative  to  the  same  radiation  level.  In  each 
case,  there  is  very  good  agreement  between  calculated  and  measured  pat¬ 
terns. 

The  above  patterns  are  presented  to  illustrate  that  our  numerical 
solutions  can  solve  the  proposed  type  of  airborne  antenna  problem.  It 
is  the  first  step  in  the  complete  verification  process.  The  next  step 
is  to  apply  the  code  to  numerous  practical  airborne  antenna  problems 
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and  ascertain  its  limits,  shortcomings,  and  possible  improvements.  As 
stated  earlier,  this  step  is  presently  being  accomplished  under  other 
contracts.  The  significant  results  of  those  findings  will  be  reported 
later.  In  any  event,  it  appears  that  the  present  solution  does  adequately 
simulate  the  engine  housings,  wing  tanks,  "T"  tails,  etc.  for  whi^h 
it  was  intended. 


Note  that  all  the  input  data  used  for  the  test  geometries  are  given 
in  the  Appendix.  In  addition,  all  the  previous  polar  plots  are  illustrated 
in  terms  of  decibels  with  10  dB  per  major  division. 
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Fiqure  14.  Comparison  of  measured  and  calculated  near-field 
pattern  for  the  test  geometry  shown. 


Figure  16.  Comparison  of  measured  and  calculated  near-field 
pattern  for  the  test  geometry  shown. 


RELATIVE  POWER  (d8)  ^  RELATIVE  POWER  <d6) 


ANGLE  (DEGREES) 


Figure  18.  Comparison  of  measured  and  calculated  near-field 
patterns  for  the  test  geometry.  Note  that  the 
engine  is  simulated  by  flat  plates. 
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APPENDIX  A 


])  Input  data  for  the  test  geometry  as  shown  in  Figure  1. 


h  ,  t- ,  r  ,  t ,  i 
'hH  .  ,  S  . 

V  3  9  lr  1  •  ,  t  ^  * 

0  •  ,  lr  '  ■ 

>lv  i  I  *|  ')  tv  1 

•  *-.  _*♦#.»«*  «  •  /-  -A  «  -A 

VI  j.  ■,  I 

0  •  ,  l  • 

1 

PI./*  lr  1*0.  I  i  A f !(  111  Alt  PLATH  A1TACMHP  TO  CYMr:OIT; 

♦  i  t  o  *  ^  O  »  f  ^  • 

I  //O/f  ,  !  •  7o  7  M  f  ~Mt  > . 

1  m’tcn  ..v  7.v> f-vv. 

I  •  /  O  /T  ‘  f  •  2  O  ( C  *  <s  ’  ^  , 

I  .7o/t  ,  I  .  /c*7Kf 

I 

I  f  A  f  1  t  T 


2)  Input  data  for  the  test  geometry  as  shown  in  Figure  4. 


ffF  ft’  t1*  d 

boy.  Mv. 

2.bf2.b,2.'jf2.b 

Id  •  .  V-  •  .  0  • 

Id  .  ,  t  ’. 

.  <’b  , .  t  « l'.  ,  ,2b  ,  j 

VO.  ,0,>o0,  I 

O . ,  (d . 

PI. A  ll:  NO.  I  ;.<l:0!Afl(  Ul.Ak  PLATH  AlTACHI:!'  TO  CYLINDHR 

A  ,i  ,0.,0.  tt  • . 

I  .7o/t,  I  .7o7H,-wO. 

1 .  7c/f,  ,-v o. 

I  .  7(WH 

I  . 7 07 1 ,  I . 7g7P,Hv  . 

PLAiH  hl-CI  ANOMLAR  PLAT!*:  ATI  ACHLiD  70  PLA7H  #  I 

A  *  r 

I  .  7  c  / 1  ,  '  .P(.7h, 
b., 7.0, -MO. 

b .  ,  7  .  t  ,  i  .i  l . 

I . 7  o 7 1 , a. po7H, ho. 

I  ♦<.,!  ,P 

I  ,  a  ,  t’ ,  ^ 
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3)  Input  data  for  the  test  geometry  as  shown  in  Figure  9. 


F ,  F ,  h  ,  i 

24.  ,y . 

2. ‘>,2  .b,2.3,2.b 

0. ,  P.  ,0 . 

3  .  b  -  ‘ )  ,14  . 

•  • , , 3 h , 3 

vid,  s-o  ,  i  ,o,3o o,  i 

id . ,  i  . 

5 

PI- Ait  NO.  IlHK.iMT  UNO  I  !  'IhL'L-N  CYLINDER  AND  ENOINF 

4  ,  i  « 0  . ,  t ( .  ,  * 1  • 

Id  .  ,  3 . 1  ,  -4  . 

Id.  ,  1  ,  -4  . 

id  . ,  7  . 1 . ,  4  . 
id .  ,  2 . l.  ,  4 . 

PL/* ‘1 1*  NO.  2;ENOIfll:  PLAIT  l-ACIMO  CYLINDER 

4  ,  t-  ,  O . ,  0 .  ,  i  . 

id.,  /.N,*4. 

1  . 3bj  ,H  .2b  I ,  — 4. 

I .333, d. 2* I ,4. 

Id  . ,  7  . !  ■ ,  4  . 

PL  A  1 1:  NO.  3  i  EfjOINP  PLATE  PARA  LLLL  TO  IT  NO 

4 ,r ,0. ,0. ,o . 

I .333,8.20 1,-4. 

I .2N3  ,  R.VOc  ,-4. 

I  .  233  ,  V . vOr , 4 . 

1 .333.8. 2  H I ,4  . 

PL/  i  I  NO.  4 1  ENOINF  IM.A’J  I-  PACINO  A  NAY  ROM  CYLIN'NR 

4  ,  r  ,  Id  .  ,  0 .  ,14  . 

I  ,2b3  ,3  .  90<  ,-4  . 
d  •  ,  I.)  .Oh/1),  —4  . 

0. ,  10  .0873 , 4 . 

I  .  33o  ,  V .  V0<  ■ ,  4  . 

PLAit  NO.  5;ttIoin  NINO  'IIP 
4,r  ,().,F.  , l  . 

O  .  ,  I O  .ob  /•> ,  —4  . 
id  . ,  I  /  .  3 ,  —  4  . 
id  .  ,  1  /  .  0 , 4  , 
id. ,  1 0. ob 
4 

3.4. 2 

3.2.4 

4.4.2 

3.4  .2 
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4)  Input  data  for  the  test  geometry  as  shown  in  Figure  15. 


r  , r ,r  ,  i 

1  (i.  ,v . 

2  .G,2  . G , 2 

vl.  ,12.  . 

0 .  ,t  . 

.  2 1> , .  t ,  (2 .  ,  .2b,  3 
112,1  12  ,.M,  R  I 

«.,D. 

6 

PI. /'ll-  NU.  I  ?1  li-'T  WI^C  M-CIION  MET WEEN  CYLINI'HR  AND  ENGINE 
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